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One and two-dimensional quantum models: quenches and the scaling of irreversible 

entropy 
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Using the scaling relation of the ground state quantum fidelity, we propose the most generic scaling 
relations of the irreversible work (the residual energy) of a closed quantum system at absolute zero 
temperature when one of the parameters of its Hamiltonian is suddenly changed; we consider two 
extreme limits namely, the heat susceptibility limit and the thermodynamic limit. It is then argued 
that the irreversible entropy generated for a thermal quench at low enough temperatures when the 
system is initially in a Gibbs state, is likely to show a similar scaling behavior. To illustrate this 
proposition, we consider zero-temperature and thermal quenches in one and two-dimensional Dirac 
Hamiltonians where the exact estimation of the irreversible work and the irreversible entropy is 
indeed possible. Exploiting these exact results, we then establish: (i) the irreversible work at zero 
temperature indeed shows an appropriate scaling in the thermodynamic limit; (ii) the scaling of the 
irreversible work in the ID Dirac model at zero-temperature shows logarithmic corrections to the 
scaling which is a signature of a marginal situation, (iii) Furthermore, remarkably the logarithmic 
corrections do indeed appear in the scaling of the entropy generated if temperature is low enough 
while disappears for high temperatures. For the 2D model, no such logarithmic correction is found 
to appear. 

PACS numbers: 


I. INTRODUCTION 

Studying emergent thermodynamics following a 
quench of a closed quantum system initially in equilib¬ 
rium with a heat reservoir has attracted the attention of 
the scientists in recent yearsir—. These studies are also 
important from the viewpoint of non-equilibrium fluctu¬ 
ation theoremsii^i^i^ and quantum Jarzynski equality^. 
Treating a sudden quench (i.e., a change) of a parameter 
of the Hamiltonian as a thermodynamic transformation 
and decoupling the system from the bath so that the 
evolution of the system following a quench is perfectly 
unitary, a generalized second law (or Clausius inequal¬ 
ity) for an isolated quantum system has been posited 
with the entropy generated in the process being defined 
in terms of the irreversible work done in the process of 
quenching^i^i^. These studies have been generalized to 
the context of Bures metric^ and Uhlmann fidelityiS, sin¬ 
gle qubit interferrometryii, open quantum systemsi^ and 
inner friction of quantum thermodynamic processes^^. 

In the above studies, the quantum system is initially 
prepared in a mixed state in equilibrium with a heat bath 
which is then decoupled and the system is subjected to 
a sudden change (or a slow ramp dictated by a protocol 
described by a unitary operator U) of one of the param¬ 
eters. Due to the stochastic nature of the work done for 
a finite system (which is not an observable^), one makes 
resort to the characteristic cumulant-generating function. 
To elaborate, let us consider a system described by the 
Hamiltonian H{X) which is at the quantum critical point 
(QCP) at A = OM- — ; we assume that at an initial time 
t = 0 the system is kept in contact with a heat bath at 
an inverse temperature j3 with an initial value of A = Aq. 
The system is then decoupled from the bath and the pa¬ 
rameter A is tuned from the initial value Aq to the final 


value A/ = Aq + (5, so that the subsequent temporal evo¬ 
lution of the system is entirely dictated by the unitary 
operator [7. Implementing the analytical continuation 
s -T —it, one defines the characteristic function 

G(s) = Tr[t/+e"^(^^)l7e-"^(^“Vo(Ao)]; (1) 

here, po(Ao) = /Z{\q) is the density matrix 

characterizing the initial Gibbs state and Z{Xq) is the 
corresponding partition function. For an instantaneous 
quench from Aq to A/, G(s) takes a simpler form: 

G(s) = Tr[e"-f^(^^)e-"^(^“Vo(Ao)]. (2) 

At this juncture, it should be noted that the charac¬ 
teristic function G(s) closely resembles the mixed state 
Loschmidt echo^^ which has been extensively studied in 
the context of zero-temperature non-equibrium dynam¬ 
ics of quantum critical systemsi^^— and associated dy¬ 
namical phase transitions^. Moreover, G(s) is related 
(through the inverse Fourier or Laplace transform) to 
the distribution function P(W) of work {W), which is 
usually characterized by an edge singularity at zero- 
temperatureSSr— . 

Focussing on the sudden quenching case, the average 
work performed is defined through the first cumulant of 
the characteristic function, as shown: 

(IF) = dG{s)/ds\s=Q = Tr {H{Xf)po) - Tr {H{Xo)po) ■ 

(3) 

To arrive at the generalized second law, we note that 
for a non-equilibrium process, the average work (IF) al¬ 
ways exceeds the free energy difference (AF) between 
the initial and final equilibrium states (both with inverse 
temperature /?) with parameters Aq and A/, respectively. 
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One can define the irreversible work through the differ¬ 
ence = {W) — AF. The heat exchange with the bath 
for a closed quantum system during the evolution being 
zero, the irreversible entropy generated in the process 
satisfies the relation AS'irr = In the subsequent 

discussion, we shall drop the superscript T and use Wirr 
to denote the finite temperature irreversible work. 

For a sudden quench of magnitude S at zero- 
temperature, when AF —>• {Eg{Xf) — Eg{Xo)), Eq. ([3]) 
reduces to 

= (^g(Ao)|i?(A/)|V;,(Ao)) - EgiXf); (4) 

here, |'i/'g(Ao)) is the ground state wave function 
corresponding to the initial Hamiltonian H(Xo) and 
Eg{Xi){Eg{Xf)) is the ground state energy of H{Xi) 
{H{Xf)). Clearly, is the same as the residual en¬ 

ergy (excess energy above the ground state) and hence 
should follow an identical scaling relation for a sudden 
quench obtained using an adiabatic perturbation theory 
as reported in the literaturei^i^Zr— . It is noteworthy that 
very recently, the quantity has been studied estab¬ 

lishing a connection to the heat susceptibility (for (5 —>■ 0) 
in the context of quantum latency which detects the or¬ 
der of a quantum phase transition (QPT)^. 

Let us now ask the question how does the proximity 
to a QCP influence the behavior of AS'irr and explore its 
scaling behavior in terms of the deviation from the QCP 
(A), the magnitude of the quench 5 = A/ — Aq and the 
system size L. In Ref. Q, it was shown how the form 
of G{s) given in Eq. ([T|) can be exploited to establish the 
Tasaki-Crooks relation and the Jarzynski equality for a 
closed quantum system. Taking the example of a ID 
quantum Ising model, it has also been shown that when 
A is varied keeping 5 small and L finite, ASirr shows a 
sharp peak at the QCP if (3 is large, while on the contrary 
in the high temperature situation (/3 —>• 0), ASirr gets 
broadened. One therefore concludes that the peak in 
ASirr is an indicator of a QPT even at a low but finite 
temperature. 

In this paper, we shall first consider and inspired 

by the scaling of the quantum fidelity explore its scal¬ 
ing behavior for a sudden quench of magnitude 6 in the 
vicinity of a QCP characterized by the associated quan¬ 
tum critical exponents in different limits A, 6 and L; both 
the situations when S can be treated as a perturbation^^ 
or it is a scaling variable will be probed. We shall then 
address the question whether ASirr generated for such a 
quench exhibit similar scaling relations. To the best of 
our knowledge, our work is the first attempt that gen¬ 
eralizes the derivation of ASirr to a higher dimensional 
system where an exact analytical calculation is possible. 
The exact analytical form enables us to derive the scaling 
of and ASirr when (5 —>• 0 and also for finite 5. 

The paper is organized in the following manner: the 
scaling relation of and AS^r are proposed in Sec. 

im In Sec. imi we shall invoke upon the one dimensional 
(ID) and two dimensional (2D) Dirac Hamiltonians to 


illustrate the scaling relations proposed earlier. The ad¬ 
vantage of using a Dirac Hamiltonian is the inherent 2x2 
nature of the same rendering it integrable in all dimen¬ 
sions; in addition, ID Dirac Hamiltonian, as we shall 
show below, provides an ideal example of a marginal sit¬ 
uation where we find the signature of logarithmic cor¬ 
rections to the scaling of which persist in AS'irr 

at low but finite temperatures and disappears in the 2D 
case. In the 2D case, in contrast to the ID case, ASirr 
does not show a sharp peak at the QCP and hence can 
not be treated as an ideal indicator of a QPT. Concluding 
comments are presented in Sec. IlYl 


II. SCALING RELATIONS OF Wirr AND ASirr 

Before embarking upon the study of ASirr at finite 
temperature, we piece together the scaling of the defect 
density and the residual energy at zero temperature fol¬ 
lowing a sudden quench considering different limits: we 
first consider the situation when L is the largest length 
scale of the problem which implies L ^ S~'', A”'^, where 
v is the correlation length exponent associated with the 
QCP; clearly, in this case the parameter 6 is finite and 
hence can not be treated as a perturbation. We shall 
arrive at the scaling relation of using heuristic ar¬ 

guments inspired by the scaling of the ground state quan¬ 
tum fidelity^i-i^, E(A, <5) = |('0g(Ao)|^/’g(Ao -f 5))| given 
byM 

\nF =-L^\S\’'^F . (5) 

where F is the corresponding scaling function, v is the 
correlation length exponent associated with the QCP and 
d is the spatial dimension; L being the largest length scale 
of the problem, this limit is known as the thermodynamic 
limit. Let us first consider the situation, A 5, which 
means 5~'' <C X~'^, i.e., the correlation length ^ ~ X~'^ is 
larger than the length scale associated with the parame¬ 
ter 5, so effectively the system is close to the QCP. Using 
the finite size scaling argument^, we expect In E to scale 
with S. Demanding that right hand side of Eq. m is 
dimensionless, we find a characteristic momentum scale 
k ^ 5'^ (since we know that has the dimension of 
length). Using the phase space argument, one then finds 
that the defect (or quasi-particle) density (n) scales as 
k‘^ ^ (For V = d = 1, it has been proved^ that 

one has n ^ |(5|.) We note that as predicted in Ref. 
[13, the scaling of defect density is different from the 
residual energy for a sudden quench in the vicinity of 
the QCP; in this case, each excitation carries energy 5'^^^ 
where z is the dynamical exponent, resulting in the scal¬ 
ing ~ 

We now consider the other limit, when X ^ 6, (imply¬ 
ing X~‘' <C S~‘', i.e., away from the QCP), A is expected 
to play the role of the scaling variable. The quantum 
fidelity defined in Eq. (O is expected to scale as InF = 
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—in a similar spirit identifying the charac¬ 
teristic momentum scale, it is straightforward to arrive 
at the scaling of the defect density as n ~ ~ 

which has already been established for the quantum Ising 
case^ {v = d = z = 1) enabling us to arrive at the scal¬ 
ing relation, ~ ^ 2 yud+i/z -2 _ p\-a q, jg 

the corresponding “specific heat” exponent. 

We have illustrated above how the scaling of fidelity 
can lead to the scaling of in the thermodynamic 

limit. Concerning the scaling of fidelity, it is well-known 
that one can switch from the thermodynamic limit to 
the fidelity susceptibility (xf(A)) limi t (where the 
notion of the fidelity susceptibility is meaningful) when 
(5“^' ^ L. Continuing along the same line of arguments, 
one can define the heat susceptibility xe^ such that 
^ S'^Xe^ when (5““" (i.e, (5 —>■ 0) is the largest length 
scale. In this limit, one expects a scaling ^ ^^A““ 

away from the QCP and on the other hand, close to the 
QCP, L plays the role of the scaling variable in lieu of A. 
To summarize, we have: 

Heat Susceptibility limit: 5~'' is the largest length 
scale and (5(—>■ 0) can be treated as a perturbation 


scaling relation for such a thermal quench. Using a per¬ 
turbative expansion of the free energy (corresponding to 
the final parameter) valid in the limit 5 —>■ 0, a scaling 
AS'irr ~ (5^A““, has been proposed^; obviously, a simi¬ 
lar expansion is not possible in other limits mentioned 
in Eqs. m and ©. In the finite temperature case, the 
system is initially in a mixed state. Nevertheless, for 
low enough temperatures the ground state of the Hamil¬ 
tonian is maximally populated and one would therefore 
expect AS'irr to satisfy a similar scaling as on 

the contrary, these scaling relations should disappear at 
higher temperatures when the initial state deviates sig¬ 
nificantly from the ground state. Our aim here is to show 
that this indeed is the case in all the limits not only for 
<5—^0. Furthermore, we shall illustrate how to arrive at 
these scaling relations in both the susceptibility limit and 
the thermodynamic limit for an integrable Hamiltonian. 
Finally, the question whether at a low temperature AS'irr 
always shows a sharp peak at the QCP will be addressed. 

III. SCALING IN DIRAC HAMILTONIANS 
AND ROLE OF MARGINALITY 


_ ^2^2/u-{d+z) ^ ^ ^ 

The crossover between the limits close and away from the 
QCP occurs when L ~ A”'". 

Thermodynamic limit: L is the largest length scale: 


Wl=^lL ~ LyS-’^> 

~ jUd+o L>\-'^ > 5-r (7) 

Similarly here the crossover from one scaling to the other 
occurs when A ~ <5. As explained above the crossover 
from the thermodynamic limit to the susceptibility limit 
occurs when ~ L. At this point, it would be instruc¬ 
tive to note that what we present above happens to be 
the most general scaling relations of . If the sys¬ 

tem is precisely at the QCP (A”*^ —oo), one arrives at 
the scaling relations (with L or 5) proposed in Ref. 
by setting A = 5 in the scaling relations and switching A 
and 5 in the corresponding conditions in Fqs. (0). 

We would like to emphasize here that when the com¬ 
bination v{d + z) exceeds 2, the scaling of lUirr is non- 
universal. In the marginal case when v{d + z) = 2, one 
encounters additional logarithmic corrections to the scal¬ 
ing and one needs to introduce an upper cut-off in the 
momentum scale fc^ax as we shall illustrate below for the 
ID Dirac model. 

We shall now proceed to analyze the scaling of AS'irr 
generated following a sudden change of A by an amount 
6 while the system was initially in thermal equilibrium 
with a heat bath at a temperature T. Question we raise 
is that wether ASirr = /5Wirr, is likely to show a similar 


In this section, we shall derive the scaling relations 
of and address the issues raised in the previous 

section concerning the scaling of ASirr using the exam¬ 
ple of ID and 2D Dirac Hamiltonians which are very 
important from theoretical point of view in quantum 
condensed matter systems. For example, we take a 2D 


massive Dirac Hamiltonian Hjj 


fdx $t(x) 


maz - 


ifiVF{dxdx+dydy) 'I'(x), where dt(x) is a two-component 


spinor field operator describing, the effective low-energy 
degrees of freedom for electrons on a honeycomb lat¬ 
tice with unequal sublattice potentials around a single 
Dirac point (i.e., single valley graphene Hamiltonian) 
with Fermi velocity vp (set equal to unity below )^. Vari¬ 
ations of this model are ubiquitous in the field of topolog¬ 
ical insulators^i^. The edge states in a 2D topological 
insulator are described by an effective ID Dirac Hamilto¬ 
nian, whereas its bulk states are described by a 2D Dirac 
Hamiltonian and the QCP separating the gapped to gap¬ 
less phases is a 2D massless Dirac point. These models 
have turned out to be immensely useful in the studies of 
the Kibble-Zurek scaling^, sudden quenche a^^A^ , fidelity 
susceptibility and thermodynamic fidelity^, Loschmidt 
echo^ and periodic steady state reached through a si¬ 
nusoidal variation of the mass term^. Considering the 
massive 2D Dirac Hamiltonian one can rescale the units 
appropriately to obtain the 2x2 Hamiltonian describing 
the system close to a single valley as. 




m 


hx iky 


iky 


( 8 ) 


here, m is the Dirac mass and the vector k is mea¬ 
sured with respect to the corners of the Brillouin zone. 
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The scaling relations obtained for the 1 D massive Dirac 
Hamiltonian: 

-m)’ (9) 

are to be compared with the corresponding 2D case. Both 
the ID and 2D models exhibit a QCP where the mass 
term m vanishes and the dispersion at the QCP is lin¬ 
ear. These QCPs are then characterized by the critical 
exponents v = z = 1. In the sense of universality Huj 
is identical to a ID quantum Ising chain; the later can 
be recast to the form given in Eq. ® with m being the 
deviation of the transverse field from its quantum critical 
value. The ID case is conspicuous in the sense that we 
have i'{d + z) = 2. As a result, one expects a logarithmic 
correction to the scaling of remarkably this 

logarithmic scaling persists in AS'in- even at low temper¬ 
ature. On the contrary, no such logarithmic singularity 


is observed in the 2D case. 

Using Hamiltonians ([8]) and ®, we shall consider a 
sudden quench of the mass term m (which will still be 
denoted by the parameter A for consistency) from an ini¬ 
tial value A to a final value A -I- 5 and calculate the exact 
analytical expressions for and AS'irr; this would in 

turn lead to the scaling relations of these quantities both 
in the susceptibility as well as the thermodynamic limit. 

Using the eigenstates of 1-D Dirac Hamiltonian given 
by 


(Ao - V 1 / 

( 10 ) 

with corresponding eigen-energies £’(Ao)i_d = 

±-\/fc2 _|_ o, the expression for W-^^. ° 

can be derived in a closed form 


|'0fc^(Ao))lD = 


1 _ _ Oqr I Oqr 

'^irrlZ)° = 2 + (A -b 5 ^ — fcmax\/fcmax + ^ ^ 


. 27r, ,2tt 


- M\n^r + ix+sr + {x+ 6 riog 


+ A(A + 25)log 


¥ + \/(¥)Uv 


fcmax + -\/fcmax + (A -|- (5)2 


fcmax + Vfcmax + (A -b 6 ^ 

I- 


( 11 ) 


Here, we have converted the sum over k to the integral 
using a periodic boundary condition with the step size 
27r/L; fcmax plays the role of the upper cut off in the 
momentum scale, i.e., fcmix which is the shortest length 
scale of the problem (i.e., the lattice spacing). This arises 
because of the fact that the integral in Eq. dl) can not 
be extended to fc —> oo and to avoid divergence one needs 
to introduce a cut-off. Furthermore, as is evident in the 
expression in Eq. (fTTll . the presence of fcmax in fact ren¬ 
ders the argument of logarithm dimensionless. Similarly, 
one can derive an exact analytical form of the in 

the 2D case also. In Fig. o, we show the variation of 
as a function of A systematically in different lim¬ 
its. We note that when plotted on the same scale, there 
is a significant peak at the QCP in the thermodynamic 
limit while the peak almost disappears in the suscepti¬ 
bility limit. 

Using the exact form of as given in Eq. m , one 

can now derive the scaling behavior of the same in dif¬ 
ferent limits: in the susceptibility limit is the largest 
length scale that never appears in the scaling. In this 
limit, we then find: 


w^rVL ~ s^\ogi2kma./xy,s-' > L > A-" > fc-L 

( 12 ) 

~ S'^ \og{Lkmax/‘ 2 Tr);S~’' > X~’' > L> ky^^^ 

(13) 

On the other hand, in the thermodynamic limit where 
L is the largest length scale, using Eq. (ED, we get 

wI=VL - 52log(2fcmax/A);L > r" > A-" > fc'L 

(14) 

-- (5^ log(2fcmax/l5); L > A”*" > 5~'' > k~l^ 

(15) 

The scaling given in Eq. ED might look identical to 
that given in Eq. ED; however, we would like to empha¬ 
size that these two relations are derived from the exact 
expression in the different limits. These scaling relations 
in the ID and 2D cases are illustrated in Figs. and 
([21) choosing the appropriate values of the parameters in 




















5 



FIG. 1; The plot of WiTr in different limits showing peak 
at the QCP (A = 0) for a ID Dirac Hamiltonian There is 
no visible peak in the susceptibility limit [L = 100, bottom 
most) and there is a sharp peak in the thermodynamic limit 
{L = 10000, topmost). In the intermediate case, L = 2000, 
there is a significant peak. Here, 5 = 0.001 and fcmax = Stt in 
all the cases. 



FIG. 3: lFi?r°/((5^L‘^) for ID and 2D Dirac Hamiltonian 

verifies the scaling proposed in the main text, showing a loga¬ 
rithmic scaling in the marginal situation (i.e. ID Dirac model) 
which is absent in the 2D case. The main figure clearly con¬ 
firms the scaling with 5 setting A = 0.01. Inset on the other 
hand shows the expected scaling away from the QCP in the 
thermodynamic limit with h = 0.01. The parameters for both 
the figures (main and inset) are L = 1000 and fcmax = tt. 



FIG. 2: The logarithmic dependence with A (away from the 
QGP) and L (close to QCP) for ID Dirac model can be seen 
which is clearly not present in the 2D case in the snsceptibility 
limit. The main figure verifies the logarithmic scaling with L 
when A = 0.01. Inset shows the logarithmic dependence with 
A for L = 100 in ID case and no such dependence for 2D case 
with Lx = Ly = 100. In the main figure as well as in the inset 
h = 0.001 while fcmax = tt. 


different limits. We find that the leading scaling behav¬ 
ior of is given by 5^; in the ID case, there exist 

additional logarithmic corrections to the scaling with a 
scaling variable appropriate in the limit under consider¬ 
ation. No such logarithmic correction is seen for the 2D 
case which can be established using a similar exact ana¬ 
lytical form of . Remarkably as we claimed earlier, 

in Fig. the parameter 5 indeed appears as a scaling 
variable in the thermodynamic limit when the system is 
close to the QCP; this is one of our most important find¬ 


ings which was never reported before. 

Let us proceed to the finite temperature situation, and 
ask the question whether similar logarithmic corrections 
are present in A^irr. The inherent 2x2 nature of the 
Dirac Hamiltonians enables us to calculate an exact ex¬ 
pression for AS'irr also. For a single k mode, one finds 




X tanh 


/pEk{Xo)\ 






2 log 


cosm3Ek{Xf)) 

cosh{l3Ek{Xo))' 


(16) 


Integrating the expression given in Eq. m over all the 
momenta modes, one can derive the exact analytical 
forms of A^irr both in ID and 2D models. Analyzing 
those exact expressions, in the appropriate range of pa¬ 
rameter values we can extract the scaling behavior of 
A^irr which we have illustrated in Figs. o and ©■ 
What is remarkable is that even though the leading be¬ 
havior of the scaling of AS'irr is still given by for both 
ID and 2D Dirac models, the logarithmic correction to 
the scaling persists in the ID case for relatively low tem¬ 
perature while it approaches the 2D case (when there 
is no logarithmic correction even at T = 0) as thermal 
fluctuations increase. This establishes our claim that the 
proposed scaling of indeed manifests in ASirr at 

low temperature. 

Finally, we ask the question whether ASirr does always 
exhibit a sharp peak at the QCP when plotted as a func¬ 
tion of the tuning parameter A (see Fig. ®) ; in ID case 
indeed it does which resembles the observation reported 
in Fig. (1) of the Ref. Q. However, we would like to 
draw attention to the symmetric nature of ASirr on ei¬ 
ther side of the QCP unlike the transverse Ising case. On 
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(a) A5irr/-L‘* shows a 
logarithmic scaling with A in 
ID Dirac with L = 100; the 
slope approaches the 2D case 
with increasing temperature 
where there is no logarithmic 
scaling. 



(b) The same for the 
dependence with L with 
A = 0.01. 


FIG. 4: We show that in the susceptibility limit, the loga¬ 
rithmic scaling of ASirr/F'^ persists for low-temperature in 
the ID Dirac case and its comparison with that of 2D Dirac 
where no such scaling is expected. We have taken <5 = 0.001 
and kmax ~ . 



(a) The scaling of AS’irr/F'* in 
the thermodynamic limit, 
showing clear log dependence 
in ID Dirac case which is 
clearly not observed in 2D 
Dirac case with 5 = 0.01. 



(b) Scaling of ASirr/L'^ with S 
shows a clear marginal 
behavior as a function of <5 
(close to QCP) for ID Dirac 
with A = 0.01. 


FIG. 5: For ID case there is a logarithmic correction to the 
scaling of AS'irr even in the thermodynamic limit which dis¬ 
appears at high temperatures when the behavior is identical 
to that in the 2D case. Here, L = 1000 and fcmax = tt. 


the contrary in the 2D case, there is no sharp peak even 
at sufficiently low temperature. This leads us to the con¬ 
clusion that above the marginal situation, AS'irr fails to 
be a good indicator of a QPT occurring at T = 0. 


focus has been restricted to the case when a parameter 
A of the Hamiltonian which is close to its QCP (at A = 
0) is suddenly changed by an amount (5. Based on the 
three length scales of the problem, namely X~'', S~'' and 
the system size L, we have predicted the existence of 
two scaling regions, namely the susceptibility limit and 
the thermodynamic limit. Inspired by the scaling of the 
ground state fidelity in these limits, we have proposed 
the scaling of in these limits both close to the 

QCP and away from the QCP; furthermore, a logarithmic 
correction with the parameter 5 appears in the scaling 
close to the QCP establishing that in the thermodynamic 
limit 5 obviously plays the role of a scaling variable. 

We then asked the question whether these scaling re¬ 
lations survive in the case of a sudden quench in which 
the initial state is a mixed state in thermal equilibrium 
with a heat bath and manifest in AS'irr- To address this 
question in a transparent way, we have used ID and 2D 
Dirac Hamiltonians where arriving at the exact expres¬ 
sions for both and ASirr is indeed possible. These 

exact expressions are then analyzed to derive the scaling 
relations in different limits. Remarkably, our study es¬ 
tablishes the logarithmic correction (which is a signature 
of marginality) appearing in the scaling of for ID 

indeed survives in the scaling of AS'irr for low enough 
temperatures. There is no trace of any logarithmic cor¬ 
rections in the 2D case where the sharp peak in AShr 
plotted as a function of A also get broadened thereby 
diminishing its usefulness as an ideal detector of a QPT. 




(a) ASirr showing peak at the 
QCP (A = 0) for a ID Dirac 
Hamiltonian. The peak height 
changes with the temperature 
keeping S fixed to 0.01 


(b) ASirr versus A shows peak 
at the QCP for 2D Dirac 
Hamiltonian for <5 = 0.01 and 
changing /3. 


IV. CONCLUDING COMMENTS 

In this paper, we have studied quantum and thermal 
quenches of a closed quantum system; in particular, our 


FIG. 6: The figure above shows ASirr versus A for ID and 2D 
Dirac Hamiltonians with different P and 5 = 0.01. There is a 
sharp peak at A = 0 for ID which gets broadened in the 2D 
case. Here, P increases from top to bottom. 
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